The ferroelectric transition in YMnC>3 from first principles 



Craig J. Fennie and Karin M. Rabe 
Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854-8019 

(Dated: February 2, 2008) 

We have studied the structural phase transition of multiferroic YMnC>3 from first principles. Using 
group-theoretical analysis and first-principles density functional calculations of the total energy and 
phonons, we perform a systematic study of the energy surface around the prototypic phase. We find 
a single instability at the zone-boundary which couples strongly to the polarization. This coupling is 
the mechanism that allows multiferroicity in this class of materials. Our results imply that YMn03 
is an improper ferroelectric. We suggest further experiments to clarify this point. 

PACS numbers: 77.80.Bh, 61.50.Ks, 63.20.Dj 



The hexagonal manganites, i?eMn03 for i?e=Ho-Lu 
and Y, are a class of multiferroic materials that are si- 
multaneously ferroelectric (FE) and antiferromagnetic 
(AFM)ii In particular, YMn0 3 is FE at room tem- 
perature (RT) crystallizing in space group P63cm (Ce v , 
Z=6)£ Above «1270K it has been shown to undergo a 
transition to paraelectric (PE) P63/mmc (Dg/i, Z=2)m£ 

The sequence of phase transitions from the low temper- 
ature FE to the high temperature PE phase has been the 
subject of much debate. Based on pyroelectric measure- 
ments and lattice constants determined by x-ray diffrac- 
tion, Ismailzade and Kizhaev argued that YMn03 un- 
dergoes a FE transition from P63cm (Z—6) to an in- 
termediate nonpolar P63/mcm (Z=6) at «930K£ In 
contrast, Lukaszewicz and Karut-Kalicinska were unable 
from x-ray diffraction to resolve this intermediate non- 
polar space group^ at 1005 K and saw no clear anoma- 
lies in the lattice constants at «930K. Katsufuji et al. 
confirmed the latter results on the temperature depen- 
dence of the lattice constants and resolved the FE space 
group up to 1000K, the highest temperature they con- 
sidered^ The FE space group immediately below the 
zone-tripling transition is also supported by the struc- 
tural refinements of Lonkai et al. on the closely related 
compounds ReMnC>3, for Re = Lu, Tm, Ybi On the 
other hand, Nenert et al. recently interpreted their x-ray 
diffraction measurements to support the observation of 
the intermediate P63/mcm (Z=6) between I020K and 
I270K. 4 The presence of some kind of feature at 930K is 
suggested by small anomalies in resistivity data at w 950 
Km (also observed in HoMnOa 9 ). Thus, the nature of the 
phase right below the 1270 K transition out of the PE 
phase, a precise characterization of the behavior around 
900-1000 K, and the identification of the PE-FE transi- 
tion in this material, remain subjects of great interest. 

It has been established that to coexist with magnetism, 
the nature of the ferroelectricity cannot be that found in 
prototypical perovskite FEs such as BaTiOs— Indeed, 
first-principles calculations have shown that the usual 
indicators of a FE instability (e.g. large Born effective 
charges) are absent in YMn03 — Polarization P s = 6 
fiC/cm 2 results from the tilting of the Mn-centered oxy- 
gen octahedra and buckling of the Y-0 planes, with no 
significant off-centering of the Mn cations as would be 



characteristic of a perovskite FE— As such, YMn03 has 
been referred to as a "geometric ferroelectric."— In this 
paper, we combine a group theoretical analysis of the ex- 
perimental FE structure with first-principles calculations 
to describe the mechanism by which this "geometric fer- 
roelectricity" originates. Computation of the coupling 
between a zone-boundary instability and a soft, but sta- 
ble, zone-center phonon shows that it induces a low tem- 
perature P s 6.5 [iCjcrr? consistent with experiments. 
Identifying this zone-boundary mode as the primary or- 
der parameter (OP), we argue that there is an improper 
FE transition in YMn03 at the observed zone-tripling 
transition, w 1270K. We suggest possible origins for the 
feature that has been observed around 1000K. 

First-principles density-functional calculations using 
projector augmented-wave potentials were performed 
within the LSDA+U (local spin-density approximation 
plus Hubbard U) methods* as implemented in the Vienna 
ah initio Simulation Package (VASP). 14-15 The wavefunc- 
tions were expanded in plane waves up to a kinetic en- 
ergy cutoff of 500 eV. Integrals over the Brillouin zone 
were approximated by sums on a 4 x 4 x 2 T-centered 
fc-point mesh. Values of on-site Coulomb interaction U 
= 8eV and exchange parameter J=0.88eV were used for 
the Mn d-orbital as calculated by Ref. 16> Polariza- 
tion was calculated using the modern theory of polar- 
ization— as implemented in VASP. All results presented 
are within LSDA+U with A-type AFM order— Calcula- 
tions of phonons, structural energetics, and polarization 
were also performed within LSDA, and with both frus- 
trated AFMi& and FM order; the results are essentially 
unchanged. 

We performed full optimization of the lattice param- 
eters and internal coordinates in the reported FE space 
group P63cm. As can be seen in Table H] excellent agree- 
ment with the experimental RT structure (lattice con- 
stants within paO.5% and bond lengths within «1% of 
the experimental values) is obtained. 

Group theoretical analysis (e.g. Ref . ITsffnif) shows that 
there are three distinct symmetry-allowed phase transi- 
tion sequences connecting the high-temperature proto- 
typic phase, P63/mmc, to the low-temperature FE phase, 
P63CIT1, as shown in Fig. ^ There is a clear experimen- 
tal consensus that the zone-tripling transition occurs be- 
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TABLE I: Crystal structure of ferroelectric YMn03, Space 
Group: P63C771, Exp. room temperature: a= 6.139A, c= 
11.41A. LSDA+U: a= 6.099A, c= 11.421 
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FIG. 1: Group-subgroup sequence of allowed phase tran- 
sitions connecting the high-temperature prototypic phase, 
P63/mmc, to the low-temperature ferroelectric phase, P63cm. 
Each sequence is labeled by an irrep of the prototypic phase 
where structural distortions associated with these irreps de- 
termine possible intermediate phases. 



fore or at the FE transition, thereby excluding P63inc 
as an intermediate phase and ruling out path (1). In 
principle, path (2) could account for both the observed 
zone-tripling transition at 1270K and a possible second 
transition at wlOOOK. This transition sequence would 
be a first-order transition to the non-polar intermediate 
phase (by the freezing- in of a K\ phonon), followed by a 
second-order, proper FE transition to the RT phase (due 
to a softening of an infrared-active phonon). Path (3) 
describes a second-order transition out of the prototypic 
phase, by the freezing-in of a K3 phonon, directly to the 
FE space group. Unlike what has been proposed in the 
literature^, Path (3) does not require a second transition 
to account for FE order (as discussed below) m& 

Independent of the pathjSI the distortion that relates 
the prototypic phase to the low temperature FE phase 
can be decomposed into the symmetry-adapted modes of 
the prototypic phase as follows: (zone-center polar 
mode), K\ and if 3 (zone-boundary modes at 
and r+ (the identical representation) ^StSi This decom- 
position is shown in Table ^ for the experimental RT 
structure. The excellent agreement between theory and 
experiment for the structure, Table P results in similar 
agreement for the mode decomposition. 

The relative strength of each mode in the FE structure 
may be quantified by considering its amplitude^ Q, also 
given in Table[H] We find, as did previous first-principles 
calculations and structural analysisjii that the FE struc- 
ture is dominated by the K3 mode, with Y and O eq atoms 
undergoing the largest displacements. The relative sizes 
of the symmetry-adapted mode amplitudes in the RT FE 
phase offer a valuable clue to which of the three phase- 
transition paths is taken by YMnC>3 . In soft-mode driven 
structural transitions one expects a mode that stabilizes 
an intermediate phase at higher temperatures to be at 
least comparable in strength (not necessarily the largest) 
to any other modes appearing in the lower temperature 
phase. If Qkx were the primary structural distortion out 
of the prototypic phase, it should be retained through 
the second transition into the low temperature FE phase, 
which is not the case in YM11O3. Thus path (2) tends to 
be ruled out by the smallness of Qk ± - Finally, we con- 
sider the remaining possibility, path (3). The fact that 



Qk 3 is one to two orders of magnitude larger than Q r - 
and Qk x respectively, suggests that Qk 3 is in fact the 
primary structural distortion out of the prototypic phase 
(e.g Ref.|2l|) and thus that YM11O3 follows path (3). 

First principles (T=0) investigation of the phonons 
and the structural energetics around the prototypic phase 
also shows path (3) to be the most likely scenario. For 
the structural parameters of the prototypic phase, we 
used the minimum energy structure within P63/mmc: 
a=6.168A, c=11.23l and u 0ap =0.9166. This corre- 
sponds to a nearly homogeneous compression of the lat- 
tice constants by about 1.4% compared with those mea- 
sured by Lukaswewicz at 1280K, an underestimate typi- 
cal of LSDA calculations. Phonons at the T and if-point 
in the Brillouin zone of the prototypic phase were com- 
puted by constructing the relevant block of the dynami- 
cal matrix from Hellmann-Feynman forces. This allowed 
us to greatly reduce the number of calculations and pro- 
vided higher numerical accuracy in the determination of 
the force-constant eigenvectors. At the T-point we find 
that the three infrared-active phonons are stable with 
the lowest frequency T^TOl): to — 90 cm -1 . At the K- 
point, we computed the three K3 phonons by a similar 



TABLE II: Decomposition of Exp. Ref. Q atomic displace- 
ments (A) into symmetry-adapted modes of prototypic 
PQz/mmc. Mode amplitudes, Q, are shown in A. 
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FIG. 2: Energy as a function of (top) Q r - and (bottom) 

1 2 

Qk 3 ■ 

method and found one highly unstable mode Kail): lu 
— i 153 cm -1 . In contrast, both K\ phonons are quite 
stable with Ki(l): w = 264 cm" 1 and K\{2): uj = 522 
cm -1 . Since lattice instabilities are known in some cases 
to be sensitive to volume, we repeated these calculations 
at the experimental lattice constants; the results are qual- 
itatively unchanged. 

The stability of the and K\ modes is additional 
evidence to rule out paths (1) and (2), respectively. Ad- 
mittedly, as these are T=0 calculations there is still a 
possibility that an intermediate phase could be stabilized 
by an anomalous entropic contribution to the free energy 
at finite temperatures^ so that, e.g. P63/mcm would 
be more stable than P63cm for a range of temperatures 
below 1270K. However, in light of the strength of the 
K3 instability and as we will show, the associated energy 
scale, the transition out of the PE phase is most naturally 
driven by the K3 mode, corresponding to path (3). 

We can understand how the -^3(1) zone-boundary in- 
stability produces a polarization by calculating the en- 
ergy surface around the prototypic phase as a function 
of the relevant modes. To simplify the analysis, the fact 
that all and Ki phonons are stable in the prototypic 
phase and contribute very little to the total structural 
distortion suggests that the coupling of these modes to 
^3(1) will be small and can be ignored. Within this 
reduced subspace we expand the energy of the crystal 
including all symmetry allowed terms to fourth order in 
Qk 3 and Q T - to obtain: 

F{Qk 3 ,Q T -) = a 20Q 2 K 3 + "02<3p- + P4qQ% 3 + A)4<2p- 
+P 3 lQ 3 K 3 Q r - +P22Q 2 K 3 Ql- (1) 

Next we performed a series of first-principles calcula- 
tions where we freeze-in various amplitudes of Qk 3 and 
Q r - and compute the total energy. The energy expan- 
sion, Eq.^ was then fit using only six data points, yield- 
ing excellent agreement with the remaining ~20 data 
points. As we will show, the model energy fits the calcu- 
lated data well over the entire energy surface, justifying 
our truncation of the Taylor expansion to fourth order. 



FIG. 3: Energy as a function of Q r - at fixed Qk 3 . 



In Fig. [21 we show the total energy (of the 30-atom cell) 
as a function of the Q r ~ (top) and Qk 3 (bottom) frac- 
tional displacements. As expected, the mode is stable 
while the double-well potential of the K3 mode is appar- 
ent. The coupling between the modes can be seen by 
plotting the energy as a function of Q r ~ at fixed Qk 3 
as m Fig. Notice that once Qk 3 becomes non-zero 
the equilibrium position of Q r - (indicated by the arrow) 
shifts to a positive value and continues to increase with 
increasing Qk 3 - Also as Qk 3 increases, the curvature of 
the energy increases due to the ^22 =0.155 eV term which 

renormalizes the quadratic coefficient of Q r - . The effect 
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of the Kz mode is not to decrease the stability of the 
mode, but rather to "push" the mode to a non-zero 
equilibrium position, playing a role analogous to that of 
a field. This "geometric field" is the mechanism for fer- 
roelectricity in the hexagonal manganitesi-i-i Further, the 
strong /?3i =-0.334 eV coupling between the ^3(1) and 
r^"(T01) lattice modes, i.e. phonons, is responsible for 
the comparatively large value of P s w 6/iC/cm 2 . Indeed, 
if this coupling were to be turned off, the polarization 
generated by ^3(1) (coupling only to the electronic den- 
sity) is reduced to P s ~ 0.1fiC/cm 2 ^ 

Now let us return to the discussion of the PE-to-FE 
path. If we now identify ^3(1) as the primary OP and 
r^"(T01) as a secondary OP, inspection of Fig. [3] reveals 
interesting crossover behavior of the mode. Minimiz- 
ing the free energy over Q r ~ leads to two cases: 

Region 1: Q Ks « y^gf => Q r _ ex -fe<&, 

Region 2: Q Ks » => Q r _ ex ~^Qk 3 

Since the polarization, P, is a first order polar tensor it 
has the same transformation properties as the symmetry- 
adapted mode, i.e. P oc Q r - . Using the modern 
theory of polarization we can calculate the polarization 
as a function of the structural mode. The physical 
quantity of interest though is the polarization as a func- 
tion of the primary OP, Qk 3 > which we show in Fig. 0] 
As the K3 structural distortion is initiated, P s becomes 
non-zero but is small due to the Q% 3 dependence until 
the emergence of a crossover to P s oc Qk 3 ■ This leads to 
an apparent "turn-on" at a finite value of Qk 3 ■ We calcu- 
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FIG. 4: Polarization, P s , vs. primary order parameter, Qk 3 - 

late the total energy (per formula unit) lowering from the 
PE to FE to be O(1000K) (the precise value fortuitously 
being 1240K) and the width of this crossover region to 
be O(100K). 

In this paper, we have shown that the overwhelming 
weight of the experimental and theoretical evidence sup- 
ports path (3), with an improper FE transition out of 
the high-temperature PE phase at 1270 K. The K 3 (l) 
phonon is strongly unstable and couples to the T^TOl) 
phonon leading to the observed polarization. The un- 
usual dependence of the polarization on the primary or- 
der parameter would lead to a broadened pyroelectric 
peak lower than T c and could, through coupling to strain, 



lead to an isostructural, FE-to-FE, transition, although 
this last point is only speculative. This differs from pre- 
vious models^ in that YM11O3 is already polar when this 
isostructural transition occurs. An additional intrigu- 
ing possibility to explain the apparent lower temperature 
transition is that, due to the presumably inhomogeneous 
nature of K3 near T c (which we have been assuming to 
be homogeneous) , fluctuations suppress the coupling be- 
tween ^3(1) and r^"(T01) till some temperature lower 
than T c , this second transition temperature likely being 
sample dependent^ In this scenario, the path to the FE 
phase is still along (3), but the P63cm phase in the in- 
termediate temperature range would be paraelectric in 
an average sense, as originally suggested by Lukaszewicz 
and Karut-Kalicinska^ Clearly, further theoretical and 
experimental work is necessary to characterize the lower 
temperature isostructural transition and the nature of 
the polarization in the intermediate phase. Finally, it 
should be appreciated, that regardless of the actual path 
taken, the softening of the K 3 mode sets up a "geomet- 
ric field" that induces a substantial polarization, thus 
providing a mechanism to achieve multiferroicity in the 
hexagonal manganites. 
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